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We present the results of magnetization, ESR and NMR measurements on single crystal samples of 
the frustrated S = 1/2 chain cuprate LiCu202 doped with nonmagnetic Zn'^^. As shown by the x-ray 
techniques the crystals of Li(Cui_a:Zna:)202 with x < 0.12 are single-phase, whereas for higher Zn 
concentrations the samples were polyphase. ESR spectra for all monophase samples {0 < x < 0.12) 
can be explained within the model of a planar spin structure with a uniaxial type anisotropy. The 
NMR spectra of the highly doped single crystal sample Li(Cuo.9Zno. 1)202 can be described in the 
frame of a planar spin glass like magnetic structure with short range spiral correlations in the crystal 
(a6)-planes with strongest exchange bonds. The value of magnetic moments of Cu^^ ions in this 
structure is close to value obtained for undoped crystals: (0.8 ± 0.1) hb- 

PACS numbers: 75.50.Ec, 76.60.-k, 75.10.Jm, 75.10.Pq 



I. INTRODUCTION 



LiCu202 is an example of 5* = 1/2 magnet with frus- 
trated exchange interactions. The magnetic structure 
of LiCu202 can be considered as a system of coupled 
S = 1/2 chains of magnetic Cu^"*" ions. Within the 
chains the nearest spins interact ferromagnetically and 
the next-nearest neighbors antiferromagnetically. Inter- 
est in magnetic systems with such a type of frustration 
is stimulated by theoretical predictions of unusual mag- 
netic phases for uncoupled and weakly coupled (ID and 
quasi-lD) chain models (see for example Refs. Il|l3)- It 
was shown that a new type of magnetic ordering can be 
realized for the magnets with this particular type of frus- 
tration. Such a type of magnetic order, with zero average 
magnetic moment on each magnetic ion, but with long- 
range ordered correlations of spin components of neigh- 
boring ions, is classified as spin-nematic order 4 The pos- 
sibility to find this new type of magnetic order experi- 
mentally is a very attractive task. According to RefUl, for 
the ID model with the intrachain exchange constants of 
LiCu202 a chiral long-range order in the low- field range 
and a quasi long-range ordered spin-density wave phase 
in higher applied magnetic fields H are expected. In the 
experiment the planar spiral spin structure was observed 
at T < T/v in the low-field range in LiCu202i^'^ At higher 
fields ^qH « 15 T another magnetic phase, which is not 
identified yet, was observed.— Probably, this phase is re- 
lated to a spin-density wave phase which is predicted for 
the ID model. A spin-nematic phase for the exchange 



integrals of LiCu202 is expected in even higher magnetic 
fields which are hardly accessible experimentally. The 
investigation of the zinc-doped crystals of LiCu202 was 
reported recentlj^ and new magnetic phases in this sys- 
tem were announced. This report has stimulated the 
Electron Spin Resonance (ESR) and Nuclear Magnetic 
Resonance (NMR) investigations of the magnetic struc- 
ture of Zn-doped LiCu202 single crystals reported here. 

In the first section a short overview of the literature 
data about the crystal and magnetic structures of the 
undoped samples of LiCu2 02 is given. In the second sec- 
tion of the paper the details of crystal growth and sam- 
ple characterization are described. In the next two sec- 
tions results of ESR and NMR experiments are presented. 
These experiments have shown that the main features of 
the low-frequency dynamics of the doped samples in the 
magnetically ordered phase are the same as for the un- 
doped samples. ESR spectra for the undoped and doped 
samples can be well described using a model of a planar 
spin structure with strong easy-plane anisotropy for the 
vector n normal to the spin plane. 

The NMR study of the Li(Cui_a:Zna;)202 single crystal 
with a; = 0.1 at T < T^ shows the coexistence of static 
correlations of magnetic moments corresponding to the 
magnetic structure of the undoped material and, simul- 
taneously, a static random distribution of spin directions. 
The temperature evolution of NMR spectra can be de- 
scribed by an increase of short range correlations along 
the crystallographic c-direction with decrease of the tem- 
perature. For the temperatures up to TAr/2 the NMR 



spectra can be described in the model of uncorrelated 
(ab) spin planes. 



II. CRYSTALLOGRAPHIC AND MAGNETIC 
STRUCTURES OF LICU2O2 

LiCu202 crystallizes in an orthorhombic system (space 
group Prima) with unit cell parameters a = 5.73 A, 
b = 2.86 A and c = 12.42 A,^ The unit ceh parame- 
ter a is approximately twice the unit cell parameter b. 
Consequently, the LiCu2 02 samples, as a rule, are char- 
acterized by twinning due to the formation of crystal- 
lographic domains rotated by 90° around their common 
crystallographic c-axis. 

The unit cell of the LiCu202 crystal contains four uni- 
valent nonmagnetic Cu+ cations and four divalent Cu^^" 
cations with spin S ~ 1/2. There are four crystallo- 
graphic positions of the magnetic Cu^"*" ions in the crys- 
tal unit cell of LiCu202, conventionally denoted as a, /3, 
7 and S. 

The Cu'^"'' ions in the different positions are weakly 
magnetically coupled and form four almost independent 
systems of spin chains. Such chains formed by one of 
four kinds of Cu^+ ions (e.g. the a-position) and rele- 
vant exchange interactions^ within the system are shown 
schematically in Fig. [1^. 

The positions of Cu^+ and Li+ ions in the crystal cell 
are given by copper coordinates^: (0.876; 0.75; 0.095), 
(0.376; 0.75; 0.405), (0.624; 0.25; 0.595), (0.124; 0.25; 
0.905) and by lithium coordinates: (0.376;0.75;0.068), 
(0.876;0.75;0.432),(0.124;0.25;0.568),(0.624;0.25;0.932). 
Schematic arrangement of the Cu^"*" and Li"*" ions in 
the projection of the crystal lattice onto the ac-plane is 
given in Fig. [1)3 

The transition into the magnetically ordered state oc- 
curs via two stages at T^ = 24.6 K and Tc2 = 23.2 KM 
Neutron scattering and NMR experiments have revealed 
that an incommensurate magnetic structure is realized in 
the magnetically ordered state (T < Tci)i^i^iii The wave 
vector of the incommensurate magnetic structure coin- 
cides with the chain direction (the b-axis). The magni- 
tude of the propagation vector at T < 17 K is almost 
temperature independent and is equal to 0.827x27r/6. 
The neutron scattering experiments have shown that 
the neighboring magnetic moments along the a-direction 
are antiparallel, whereas those along the c-direction are 
parallel. The investigation of the spin-wave spectra 
by inelastic neutron scattering^ shows that the incom- 
mensurate magnetic structure in LiCu202 is caused by 
a competition between the ferromagnetic exchange in- 
teraction of the nearest-neighbor magnetic ions in the 
chain Ji = —7.00 meV and the antiferromagnetic in- 
teraction of the next-nearest neighbor ions in the chain 
J2 = 3.75 meV. The antiparallel orientation of the mag- 
netic moments of Cu^"*" between neighboring chains is 
caused by the strong antiferromagnetic interaction J3 = 
3.4 meV. The coupling of the Cu^+ magnetic moments 



along the c-direction and the couplings between the mag- 
netic ions in different crystallographic positions a, /3, 
7 or (S are much weakeri^ Thus, the magnetic struc- 
ture of LiCu2 02 in the magnetically ordered phase can 
be considered as quasi-two dimensional. The quasi-two- 
dimensional character of the magnetic interactions in 
LiCu202 was also proven by resonant soft X-ray mag- 
netic scattering cxperimentsJ^'i^ 

The magnetic structure of LiCu202 at zero magnetic 
field was studied using different experimental methods, 
but so far there is no generally accepted model for it. For 
a review of proposed models see, for example, Ref. |6|. 

In magnetic fields above 3 T, the magnetic structure 
realized in LiCu2 02 is much clearer. ESR and NMR 
studies of LiCu202 in the low-temperature magnetically 
ordered phase {T < Tc2) established that a planar spiral 
magnetic structure is formed in this compoundJ^ The 
magnetic moments located at the a (/3, 7 or 5) position 
of the crystal unit cell with coordinates x, y, z (measured 
along the a, b and c-axis of the crystal, respectively) are 
defined as: 



/.t„ = /i • li(-l)^/° • cos(fc,c ■y + (j)a) + 

+M-l2(-l)"/°-sin(fc,e-2/ + 0a), (1) 

where li and I2 are the two mutually perpendicular unit 
vectors, kic is the incommensurability vector parallel to 
the chain direction (b-axis), fi is the magnetic moment 
of the Cu^"*" ion and a; is a multiple of the chain period 
a. The magetic moment per copper ion at T < 10 K 
was evaluated as /.i = 0.85 ^,3^^ The phases (/)„, 4)15 , 
(j)-y and (f>s determine the mutual orientation of the spins 
in the chains formed by the ions in the different crystal- 
lographic positions. Their values are extracted from the 
NMR dataiiii^: (pa = 0\ (pjs = (/'a+7r/2; <j)j = TT-kic*b/2; 

05 = 07 + 7I'/2. 



III. SAMPLE PREPARATION AND 
EXPERIMENTAL DETAILS 

Single crystals of Li(Cui_a;Zna;)202 with the size of 
several cubic millimeters were prepared by the "solu- 
tion in the melt" methodJ^ Samples were grown in 
air by a flux method in alundum crucibles. The mix- 
tures of analytical grade Li2C03, CuO and ZnO of 
Li2C03-4(l-x)CuO-4xZnO compositions were melted at 
1100 °C, and then solidified by cooling to 930 °C at 
a rate of 5.0 °C/hour. Next, the crucible was with- 
drawn from the furnace and placed on a massive cop- 
per plate to ensure rapid cooling to room temperature. 
Quenching from si900 °C was necessary, because the sin- 
gle phase LiCu202 single crystals decompose below this 
temperature f^ 

The samples were shaped as a flat plates with the de- 
veloped (a6)-plane. The twinning structure of the sam- 
ples was studied with optical polarization microscopyj^ In 
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FIG. 1: a. (color online) Schematic representation of the ar- 
rangement of Cu^"*" moments in LiCu202. Only one of four 
Cu^+ (q, /3, 7, 5) positions is shown. Arrows correspond to 
the spin directions below Tc2 for H || c. Ji, J2 and J3 are 
the main exchange integrals.— b. Schematic arrangement of 
the Cu^"'' and Li+ ions in the projection of the crystal lattice 
onto the ac-plane^^. Large circles indicate the Cu'^^ ions, and 
small circles - the Li""" ions. Closed and open circles repre- 
sent ions with the coordinates equal to 0.25& and 0.75b along 
the &-axis, respectively, c. Scetch of the spin configuration 
corresponding to the extreme values of dipolar field on Li site 
for two incommensurable domains with incommensurability 
vector directed upchain (solid arrows) or downchain (dashed 
arrows). All spins vectors and external field are in plane of 
figure. 



most cases, the samples were twinned with a characteris- 
tic domain size of several microns. It was possible to se- 
lect samples without twinning structure for Zn concentra- 
tions up to X = 0.1. The absence of the twinning struc- 
ture in the samples selected for the experiments was con- 
firmed by X-ray diffraction and ESR measurements i^ii^ 
We have used samples without twinning for most of the 
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FIG. 2: The concentration dependence of the crystal unit 
cell volume (top panel) and crystal unit cell parameters a,b,c 
(lower panels) of single crystals of Li(Cui_a;Zna:)202. x — 
concentration of zinc in initial charge. The solid symbols show 
the results obtained on our samples discussed here. The open 
symbols show the data taken from Ref. [3. 



experiments described below. 

The X-ray diffraction patterns were taken in a 6* — 20 
geometry using CuKa irradiation. The samples produced 
clearly distinguishable diffraction patterns from the (hc)- 
and the (ac)-plancs of the crystal. The spectra obtained 
from twinned samples had a form of a superposition of 
these spectra. The Zn doping of LiCu202 results in a 
monotonic change of the cell parameters up to a con- 
centration limit of Xc ~ 0.12. For samples with zinc 
content of initial charge x < 0.12 no diffraction patterns 
due to impurity phases were observed. Additional x-ray 
diffraction patterns were observed for the samples with 
X > 0.12, which can be ascribed to the presence of in- 
clusions of impurity phases. The patterns ascribed to 
the main phase show that the cell parameters for higher 
Zn content do not change. The zinc concentrations of 
Li(Cui_a;Zna;)202 single crystals were measured by elec- 
tron probe micro analysis on the spectrometer "Eagle 
11" ("EDAX", USA). The concentration values of x in 
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FIG. 3: (color online) Temperature dependencies of the mag- 
netic susceptibility M{T)/ ^qH for field directions H || a, b 
(lower panel) and c (upper panel) and ^.qH = 0.1 T for 
a twinned doped single crystal of Li(Cui_j;Zna;)202 with 
X = 0.1 (black lines) and for the undoped untwinned crys- 
tal a- = at H II a, b and c taken from Ref. [^. Dashed 
lines show the same M{T)/ij,oH curves for undoped samples 
(a; = 0), scaled on the y-axis by a factor of 1 — 2a; = 0.8. 
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FIG. 4: (color online) Temperature dependencies of the mag- 
netization derivative dM/dT for field directions H || a, b 
(circles) and c (squares) and ^oH = 0.1 T for twinned doped 
single crystals of Li(Cui_i,Zni)202 with a; = 0.1 (lower panel) 
and for the undoped untwinned crystal a; = (upper panel, 
taken from Ref. [y). Data are obtained from M(T)/ ^.^H at 
iioH =0.1 T. 



our single-crystalline samples obtained with this method 
coincide with the values xq of the zinc concentration in 
the initial charge with precision better than xo ± 0.03 
for all samples with x < Xc- 

The dependence of the crystal cell parameters a,bfi of 
the single crystals of Li(Cui_a;Znj^)202 on zinc concen- 
tration X in initial charge is shown in Fig. [2] The solid 
symbols show the results obtained from the samples de- 
scribed above. In contrast with the data of Ref. [3 (open 
symbols), we did not observe any anomaly in the concen- 
tration dependence of the cell constants at a; = 0.055. 

Magnetization curves in static magnetic fields of up to 
7 T were measured with a commercial SQUID magne- 
tometer (Quantum Design MPMS-XL7). 

The ESR experiments were performed with a 
transmission-type spectrometer using resonators in the 
frequency range 18 < i^ < 120 GHz. The magnetic field 
of a superconducting solenoid was varied in the range of 



< /.to-ff < 8 T. Temperatures were varied within the 
range of 1.2 < T < 30 K. 



The NMR experiments were performed with a con- 
ventional phase coherent, homemade pulse spectrometer 
at a fixed frequency of i' = 61 MHz. We investigated 
the '^Li (/ = 3/2, 7/27r = 16.5466 MHz/T) nuclei using 
spin-echo technique with a pulse sequence 5/xs-td-10//s, 
where the time between pulses, r^i, was 40 iis. The spec- 
tra were collected by sweeping the applied magnetic field 
between 3.5 < jjlqH < 3.9 T. The temperatures were 
stabilized with a precision better than 0.02 K. 



IV. EXPERIMENTAL RESULTS 



A. Magnetization measurements 



The experimental curves M(T)/ hqH and their tem- 
perature derivatives for different field directions H 1 1 a, b 
and H II c at hqH = 0.1 T are shown in Figs. |3] and 
m The dependencies were obtained on the single crystal 
without doping {x = 0) and on the twinned single crystal 
with X = 0.1. For all field directions broad maxima of 
M{T)/^qH were observed at T = 38 K for samples with- 
out doping and at T = 31 K for doped samples with x = 
0.1. The maximum is typical for low-dimensional anti- 
ferromagnets. In the paramagnetic region, the magnetic 
susceptibility for H || c exceeds that for H || a, b, which 
is consistent with the anisotropy of the g-tensor measured 
in the ESR experimentsii {ga.b = 2.0; gc = 2.2). At high 
temepatures T > 150 K the susceptibility of the doped 
sample is (20 ± 1.2)% less then that of the pure sample. 
This reduction is close to the factor (1 — 2x) expected 
if all Zn^+ ions take the places of the Cu^+ ions, which 
would result in the corresponding decrease of the num- 
ber of paramagnetic spins. To emphasize this fact the 
M{T)/^oH dependency for the sample without doping 
reduced by factor 0.8 is shown in the same figure with 
dashed lines. This experiment allows to conclude that 
the number of magnetic Cu^"*" ions in the doped samples 
of Li(Cui-xZn2;)202 is reduced by the number of non- 
magnetic Zn^+ ions. 

The transitions to the magnetically ordered states 
at Tci and Tc2 are marked by the change of the 
M{T)/^qH slope. These inflection points are well re- 
solved in the temperature derivative of M{T)/ fx^H for 
the undoped sample (Fig. 15). The {d/dT){M{T)/ ^i^H) 
curves are strongly anisotropioS: for H || b and c, 
there are two anomalies corresponding to two transi- 
tions at Tci and Tc2- The temperature dependence of 
{d/dT){M\T)/iJ.oH) for H || a shows only one sharp 
peak at the lower temperature near Tc2 in the studied 
field range. For the doped sample the inflection points 
on the M{T)/fioH curves were smoothed, but the max- 
imum of the slope for H || c was observed at evidently 
higher temperature than for H || a, b. The distance be- 
tween the maxima for doped samples with x < 0.1 are al- 
most the same as for the samples without doping (Fig.|3|)- 
Thus the double-step transition in the magnetically or- 
dered phase also takes place for doped samples. Note 
here that the double-stage transition is specific for a pla- 
nar spin structure with a strong easy-plane anisotropy for 
the vector n normal to the spin plane. The anomaly of 
M{H) at a field of /ioilci ~ 2 T for H || b which was as- 
sociated with the reorientation of the spin plane observed 
earlier for samples without doping^ was not observed for 
all studied samples with x > 0.025. 
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FIG. 5: (color online) Upper panel: ESR frequency-field de- 
pendency for Li(Cui-2,Zna;)202 samples with x — 0.075, x = 
0.1 (black and blue solid symbols), and a; = (open symbols, 
from Ref[l3), H || c, T = 4.2 K. Solid and dashed lines show 
fitting of the spectra with the equation given in the text. The 
dash-dotted line shows the spectra of the absorption compo- 
nent "C" , which coincides with the EPR spectra measured at 
T > Tjv. Lower panel: Example of ESR absorption spectra 
measured on three samples with x — 0.1, x = 0.075, and a; = 
at close frequencies {v « 79 GHz), H 1| c and T = 4.2 K. In- 
sert to the lower panel: dependence of the "A"-component 
ESR line width on the Zn concentration. H || c, T = 4.2 K, 
i/ ^ 79 GHz. 



B. Electron Spin Resonance experiments 

The ESR in single crystals of Li(Cui_a;Zni,)202 was 
studied on the samples with x — 0.025,0.05,0.075,0.1. 
The samples with x = 0.025 and a; = 0.1 were untwinned. 
ESR on the sample without doping was studied in a pre- 
vious work.— The typical ESR absorption lines are shown 
in the lower panels of Figs. [5] and ID We present the lines 
for the samples with the Zn concentration x = 0.075, 
x = 0.1 and the undoped untwinned sample from Ref. Il4l. 
The ESR absorption lines are systematically broader for 
the samples with larger Zn concentration. The insert in 
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FIG. 6: (color online) Upper panel: ESR frequency-field de- 
pendencies for Li(Cui_a;Zn2:)202 samples with x = 0.075, x — 
0.1 (black and blue solid symbols), and a; = (open symbols, 
from Ref flJI) . Intense components are marked by large sym- 
bols, weak components are marked by small symbols, respec- 
tively, for H II b, r = 4.2 K. The solid and dashed lines are 
guides to the eyes. The dash-dotted line corresponds to the 
absorption component "C", which coincides with the EPR 
spectra measured at T > Tjv. Lower panel: Example of field 
dependencies of the transmitted power, measured on three 
samples with x — 0.1, x = 0.075, and a: = at two close 
frequencies (u f» 79 GHz), H || b and T = 4.2 K. 



the lower panel of Fig. [5] shows the concentration depen- 
dence of the absorption line width for H || c at T = 4.2 K. 
The ESR spectra iy{H) for field directions H || c, b for the 
samples with x ~ 0, 0.075 and 0.1 are shown in the upper 
panels of Fig. [5] and [51 The sample with x = 0.05, 0.075 
were twinned, but due to strong anisotropy (see below) 
it was simple separate the absorption lines only from the 
domain with H || b. The dash-dotted lines in the upper 
panels of Fig.[5]and[n]depict spectra of electron paramag- 
netic resonance for LiCu202: i^ = {gnB/h)H {ga,b = 2.0; 
Qc = 2.2).— Resonance fields of the relatively weak com- 
ponent "C" are well described by such a paramagnetic 
dependence. Probably, this component is due to param- 
agnetic inclusions. 

The experimental data which correspond to the ab- 
sorption component "C" are not shown in the upper pan- 
els of Fig. Eland El The ESR spectra (i^(i?)) of the in- 
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FIG. 7: (color online) (ab)-plane angular dependencies of the 
resonance fields for "A" and "B" components measured on 
Li(Cui-a;Zna;)202 samples with x — 0.1 (solid symbols), and 
X — (open symbols) for frequencies around 80 GHz. T = 4.2 
K. The dashed lines are guides to the eyes. 



termediate Zn concentrations x = 0.025, x = 0.05 (not 
shown) and x = 0.075 interpolate between those of the 
samples with extreme doping concentration a; = and 
X = 0.1. The error bars on the Figs. 5, 6 show the antifer- 
romagnctic resonance line widths measured at the half of 
the absorbed intensities. 

The main features of the ESR spectra {v{H)) for dif- 
ferent Zn concentrations remain the same as for the un- 
doped sample. The ESR spectra {v{H)) consist of a sin- 
gle branch for H || c and two branches for H || b. The gap 
of the spectra v{H) grows from v{H = 0) = 30 ± 2 GHz 
for the undoped sample {x = 0) to v{H = 0) = 42 ± 2 
GHz for doped samples with x — 0.1. Rotation of the 
static field into the a-direction flattens the v{H) de- 
pendence, which flnally becomes fleld independent for 
H II aw^ Angular dependencies of the resonance fields 
on rotation in the (a6)-plane for samples with a; = 0, 0.1 
are shown in Fig. [71 Both resonance fields shift to higher 
fields as H approaches H || a. Such a behavior was ob- 
served also for a rotation of the applied magnetic field H 
into the (ac)-plane. From these data we conclude that 
the strong uniaxial anisotropy along the a-direction re- 
mains for all studied doped samples. Such an anisotropy 
corresponds to the easy (6c)-plane anisotropy for the vec- 
tor n normal to the spin plane. 

The branches marked as "A" for H || b,c arc quasi 
linear with field: u — k\/ H'^ + A^. The coefficient k is 
noticeably smaller than the gyromagnctic ratio gfis/h. 
Such a field dependence is typical for a planar spin struc- 
ture with strong "easy-plane" anisotropy for vector n per- 
pendicular to the spin plane and field direction perpen- 
dicular to the anisotropy axis (H _L a)J^ The coefficient 
k is deflncd by the anisotropy of the susceptibility of the 
spin structure: k = {gfJ.B/h)^y{xJJ)(J~^^T) (here x\\ and 



X± are the susceptibilities for field directions parallel and 
perpendicular to the vector n of the spiral spin structure). 
Using the value of k as a fit parameter we obtained that 
the anisotropy of the susceptibility x\\/x± ~ 1 slightly 
decreases from 0.55 ± 0.02 to 0.5 ± 0.02 with the increase 
of doping from a: = to x = 0.1. 

For the field orientation H || b the magnetic structure 
of the undoped LiCu2 02 samples undergoes a reorien- 
tation at a field of fioHd ~ 2 T. Earlier experiments 
have shown that at this field the spin plane of the spi- 
ral structure rotates from the {ab) to the (ac)-planei^ii^ 
The magnetic susceptibility at this field undergoes a step 
like increase. The ESR frequency-field dependence for 
this orientation is shown in the upper panel of Fig. |6l 
The transition field Hd is marked by the jump of the 
frequency-field dependence for the undoped sample. The 
intensive component of the ESR absorption line changes 
form the "B" -branch of the frequency-field diagram to 
the "A" -branch at the transition field Hd for the un- 
doped sample. Thus, the "A" -component dominates the 
absorption spectrum of the undoped sample at high fields 
(Figure ini). For the doped sample with x = 0.1 the ab- 
sorption component "B" was more than 10 times more in- 
tensive than the component "A" in the whole field range. 
No jumps of the position of the intense component in the 
frequency-field diagram was observed for the doped sam- 
ples. 

From these observations we conclude that for doped 
samples no magnetic reorientation takes place in the 
studied field range. This conclusion is in agreement with 
the absence of anomalies in the M{H) curves mentioned 
in the previous subsection. Note that the value of Hd 
and the sharpness of the transition are strongly depen- 
dent on the quality of the samples even for the crystals 
without dopingi^ii^ 



C. Nuclear Magnetic Resonance experiments 

The magnetic structure of the single crystals of 
Li(Cui_xZn:c)202 with x = 0.1 was tested with the NMR 
technique. ^Li NMR field-sweep spectra measured at dif- 
ferent temperatures at field orientation H || c and fixed 
frequency i/ = 61.0 MHz are shown in Fig. [5] The shape 
of the NMR spectra transforms from a single compo- 
nent spectrum at high temperatures T > Tjv to a spec- 
trum with characteristic shoulders at temperatures below 
Tjv ~ 18 K (sec Fig. 8). The spin echo NMR signal in 
the low temperature range was observed in a broad field 
range fioA.H w 0.25 T which is approximately the same 
as for LiCu202 without doping (see Fig. [3]). At temper- 
atures close to Tjv the central line of the NMR spectra 
dominates. This fact demonstrates, that by nearing to 
Neel temperature the number of nuclei with paramag- 
netic environment increases. 

In Fig. ini the lines show the fitting results of the NMR 
spectra. The ^Li NMR spectra for doped and undoped 
samples were fitted by taking into account dipolar fields 
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FIG. 8: (color online) NMR spectra of ''Li nuclei in 
Li(Cui_a;Zni)202, X = 0.1, measured at different tempera- 
tures. H II c and pumping frequency u — 61.0 MHz. 



within the model of a planar spiral magnetic structure 
in the (a6)-plane with an incommensurate wave vector 
kic = 0.827 X 27r/6 aligned along the b-direction and an 
effective magnetic moment on the Cu^"'' position equal to 
0.8/iB. Note, that the contribution of the contact field 
for such spin arrangement (H || c) can be neglected. The 
magnetic structure for the sample without doping is de- 
scribed by Eqn.(l). The mutual orientation of the spins 
from neighboring (a&)-planes defined by phases 4>a,i3,-y,s 
are given in the end of the section II. The line width 
of individual groups of resonating lithium nuclei for this 
computation was taken as 5H = 0.005 T, which is much 
smaller than all peculiarities on the experimental spec- 
tra. This value of 5H roughly corresponds to the value 
of quadrupolar splitting of NMR line on ^Li nuclei in 
LiCuaOa {vq ~ 51.7 kHz, Ref.i). 

To explain the main features of our experimentally ob- 
served NMR spectra we consider effective fields at lithium 
sites which are generated by an individual magnetic chain 
of Cu^+ ions nearest to the lithium nuclei (see the frag- 
ment of the spiral structure shown in Fig. Ic). For sim- 
plicity we suppose that the Li^ ion, the Cu^"*" chain, 
the spin vectors, and the applied static field H lies in 
the plane of the figure. The projection of the effec- 
tive field from neighboring spins on the applied field H 
changes along the spin chain sinusoidally. For a given 
propagation vector kic the NMR spectra exhibit a shape 
with two characteristic maxima at the extremum fields 
ly/j ± Hextr- Since the Li"'' ions are located symmetri- 
cally between two Cu^"'" ions, the extreme fields H^xtr will 
correspond to the arrangement of the spins in the cop- 
per chain symmetric with respect to the chosen Li ion 
(marked with index i in Fig. Ic). As it is schematically 
shown in Fig. Ic, such a symmetric arrangement corre- 
sponding to the magnetic structures, with opposite wave 
vectors kic (solid arrows) and — fc^c (dotted arrows), yields 
different dipolar contributions to effective field. Thus for 
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FIG. 9: (color online) Squares represent NMR spectra of ^Li 
nuclei in Li(Cui_a;Zna;)202, x ~ 0.1 (top and middle panels) 
and a; = (bottom panel) measured at T < Tn for H || c and 
V = 61.0 MHz. The solid lines show the simulated spectra 
in models of spiral planar spin structures within the (ab)- 
plane with long-range order in the c-direction for the undoped 
sample (bottom panel) and with disorder in the c-direction 
for the doped sample (top and middle panel). The dashed 
lines show the simulated spectra in the models of randomly 
disturbed planar spiral spin structure within afo-planes with 
short-range order in c-direction (top and middle panels). 



such simple model structure we can expect a broad NMR 
spectrum with four characteristic maxima corresponding 
to two magnetic domains with opposite wave vectors. 
The magnetic structure expected for LiCu202 at an 



arbitrary orientation of the applied magnetic field pro- 
vides NMR spectra with 8 maxima corresponding to the 
four extremal fields at the lithium nuclei from different 
positions Lil, Li2, Li3, and Li4 (Fig lb). For the field 
orientation H || c the number of maxima of the spectra 
is reduced to four (sec bottom panel of Fig. [S]). 

In order to interpret the NMR spectra of doped sam- 
ples we simulated the NMR spectra for the structures 
with random static deviations of the spin directions from 
the spin orientations of undoped model. For modeling 
purposes the deviation angles within the spin plane were 
set to the random value — 0('""^') < < 0("'°^), where 
the maximal deviation (jjmax was used as a model fitting 
parameter. Taking into account that the exchange inter- 
actions between planes is weaker than the in-plane inter- 
actions, we used different deviation parameters (j).!^""^ 

for spins within (a6)-plane, and 0^"^-' for spins from 
different planes. The NMR spectra of doped samples at 
9 K can be fitted in the model of totally random mu- 
tual orientation of the spins of the neighbor planes, i.e. 
^{max} ^ -^gpo (ggg bottom panel of Fig.9). The NMR 
spectra obtained in this model are only weakly sensitive 
to the disorder within the individual spin plane. The 
solid blue line shows the computed NMR spectra with 
the assumption of ordered moments within each (ab)- 

plane (0[™"''^ =0°). The doted green line shows NMR 
spectra obtained under the assumption that the spin di- 
rections of ions within the (a6)-planes are oriented ran- 
domly around undisturbed directions within the angle 
Arnax) _ gj-o^ p^^. ^j^^ ^^ ^^ took 0.8/ub as the absolute 

value of the magnetic Cu^^ moments at this particular 
temperature of T = 9 K. Decreasing the temperature the 
central maximum of the spectra start to broaden (see 
Fig. 8). Such transformation can be explained by the 
appearance of increasing spin correlations between ions 
of neighboring (a6)-planes towards lowering the temper- 
ature. The NMR spectrum measured at 1.8 K (squares) 
and the fitting results (lines) are given in the top panel 
of Fig. 9. Again, the solid blue line shows the computed 
NMR spectrum with the assumption of ordered moments 
within each (ab)-plane (0y°^ = 0°) and random devi- 
ations of nearest spins from neighboring planes within 



the angle 



;(max) 



45°. A better agreement between 



our model and the low-temperature experimental NMR 
spectrum is obtained for the case of a magnetic struc- 
ture with short ranged static correlations for both in- 
plane and inter-plane correlations according to the green 
dashed line. It shows the computed NMR spectrum with 
the assumption of random deviations of spins described 
by the parameters (/y^™"""^' = (j)^™""^' = 35°. Again, the 

absolute value of the magnetic moments of Cu^+ was 
taken to be 0.8/is for the fit. Here we must note that 
the local field at the probing nuclei is defined mostly by 
four nearest coordination spheres, which means that the 
NMR spectra are sensitive to short-range correlations 
only. These correlations must be static at least during 



the time window of w 0.1 s set by our NMR experiment. 



V. DISCUSSION 

X-ray diffraction and electron probe micro analysis 
of the single-phase crystals of Li(Cui_3;Zn^)202 show 
that Zn ions enter into the lattice without destruction 
of the crystal structure in the broad concentration range 
< a; < 0.12. The scaling of the high-temperature sus- 
ceptibility as (1 — 2x), shown in the Figure 3, indicates 
that Zn ions substitute for Cu^+ and each Zn^+ ion re- 
duces the number of magnetic ions by one. Nonmag- 
netic dilution should also affect the effective interspin 
interactions (i.e. Weiss temperature Q). The fit of our 
H\\c susceptibility data in the 150K — 30QK tempera- 
ture range shows that the Weiss temperature decreases 
from (61 ± 3)K for the pure compound to (51 ± 5)K for 
the 10% doped compound. The Weiss temperature of the 
doped system is also about 80% of that for the pure com- 
pound, which is again close to the (1 — 2x) factor describ- 
ing the decrease of the average number of neighbors of 
the magnetic ion with doping. On the other hand, even 
for the 10% doping the M{T) curves demonstrate the 
same characteristic broad maximum around 35i^ as for 
the pure compound. Thus, magnetization measurements 
prove that the magnetic subsystem of Li(Cui_2:Zn2:)202 
can stiU be considered as a system of spins {S = 1/2) 
with frustrated exchange interactions, where the spins 
are arranged in chains along the crystallographic b-axis. 

According to Ref. the strongest exchange interac- 
tions for undoped Li(Cui_3;Znx)202 are the interactions 
between copper ions within the (a6)-planes: intrachain 
interactions of nearest spins Ji, intrachain interactions 
of next-nearest spins J2 and interchain interaction J3 
(see Fig.l). These interactions lead to a spiral long-range 
magnetic order with an incommensurate wave vector di- 
rected along the crystallographic b-axis. The solitary 
nonmagnetic Zn^+ defects disturb this spiral magnetic 
structure in their immediate surroundings only. 

For small concentrations x the characteristic size of 
the disturbed area around nonmagnetic defects can be 
evaluated as ^ \ ^' cell units, which yields a value of 4 
cell units in the case of LiCu202. For x = 0.1 the average 
distance between such nonmagnetic defects within the 
chains amounts to « 5 lattice constants and, therefore, 
any long-range order of the magnetic structure is most 
likely destroyed for such highly doped compound. 

We suggest that such a disordered state can be a re- 
alization of a novel spin-glass like state with short-range 
chiral correlations in the (a6)-planes. It is important that 
for the system with frustration of exchange interactions 
discussed here, the replacement of magnetic ions by non- 
magnetic does not break the exchange bond within the 
chains, because such replacement excludes nearest but 
not next-nearest bonds. This circumstance can be the 
reason of the relatively weak reduction of the ordering 
temperature by doping of the samples observed experi- 



mentally. 

The ESR experiments show that the main fea- 
tures of the low-frequency spectrum of excitations in 
Li(Cui_2:Zn2-)202 do not change at a doping of a; < 0.1. 
The spectra can be explained in the frame of a planar 
spin structure with strong easy-axis anisotropy along the 
a-axis of the crystal. This means that the spins of this 
spin-glass like structure lie in the plane defined by the 
applied magnetic field H and the axial-anisotropy vec- 
tor along the a-axis for all concentrations x within the 
entire range < cc < 0.1. The absorption line width 
of Li(Cui_2;Zna;)202 strongly increases with the Zn con- 
centration (insert in Fig. 4) and the absorption lines 
of the highly doped samples are asymmetric: the low- 
field part of the absorption lines is more broadened than 
the high-field part. Probably, this asymmetry is due to 
a nonuniformity of the magnetic structure of the doped 
Li(Cui_2:Zn2^)202. The defects in the magnetic structure 
might trigger excitations of rnagnons with nonzero wave 
vectors (see for example Ref. [l^). The resonance condi- 
tion for such excitations is expected for applied magnetic 
fields lower than the antiferromagnetic resonance field. 

The NMR experiments indicate that the value of the 
magnetic moments at the position of Cu^+ ions at low 
temperature in doped samples is nearly the same as in 
the samples without doping at temperatures less then 
15 K. Additionally, the shape of the NMR spectra of 
zinc doped samples indicates short-range static correla- 
tions the same as in undoped samples. These correla- 
tions appear at the temperature « 18 K, where was ob- 
served the anomaly on magnetization curve and the ESR 
spectra start to be gaped. The temperature evolution 
of the shape of NMR spectra of doped samples can be 
explained by appearance of short-range correlations be- 
tween spins of neighbor (a&)-planes with decrease of the 
temperature from 9 to 1.8 K. The interpretation of the 
experimental spectra and our fitting results is natural 
and one of the simplest. But we can not exclude other 
models, which can describe the NMR spectra in another 
way. The results of NMR experiments on strongly doped 
Li(Cui_a;Zn^)202, for x — 0.1 are in agreement with the 
spin-glass like magnetic structure suggested above. 



VI. CONCLUSIONS 

Untwinned single-crystalline samples of 

Li(Cui_2;Zn^)202 were grown for < x < 0.12. 
It is shown that the zinc doping diminishes the number 
of magnetic copper ions Cu^"*" and a single crystallo- 
graphic phase is maintained for all samples within the 
entire doping range < a; < 0.12. The ESR spectra for 
all doping concentrations within < x < 0.12 can be 
well explained in the model of a planar spin structure 
with strong easy-plane type anisotropy for the vector 
n normal to the spin plane. The NMR spectra of 
the highly doped single crystal Li(Cuo.9Zno. 1)202 can 
be well described by planar spin-glass like magnetic 
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structure with static short-range spiral correlations. 
The value of magnetic moments of Cu^"*" ions in this 
structure is close to the value obtained for undopcd 
crystals: (0.8 ± 0.1) ^s. 



VII. ATTACHMENT 

According to the theory of exchange symmetryi^ the 
order parameter in a coplanar exchange structure is rep- 
resented by two mutually perpendicular unit spin vectors 
li and I2, which transform according to a particular irre- 
ducible representation of the crystal symmetry group. In 
a helical magnetic structure, the order parameter with 
respect to the subgroup of translations transforms ac- 
cording to the representation with an incommensurate 
wave vector. Lagrange function for this system can be 
written as 



600 



L = ^((ii+7[liH] 



-(l2+7[l2H])2) + ^(n+7[nH])2-C/„ 



where n = [I1I2], and Ua - small addition to the ex- 
change part of the energy due to relativistic effects. li 
and the susceptibility tensor components are related by 
the expressions: xi^ — Ii + -^2,X|| = 2/i. In the spiral 
magnetic structure, Ua depends only on the orientation 
of the vector n with respect to the crystal axes. In the 
simplest case, when the direction of all the spirals coin- 
cide, senior terms in the expansion of the quantity Ua 
in the components of the vector n for a biaxial crystal 
can be written as ^{Ani + Brii). For a magnetic field 
directed along one of the symmetry axes, vector n also 
directed along one of the twofold axes. The coordinate 
axes are chosen so that, in a zero magnetic field, the vec- 
tor n is parallel to the z-axis. Then, the coefficients A 
and B satisfy the inequalities A < Q, A < B. In the 
absence of the magnetic field, oscillation frequencies are 
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The third oscillation frequency is equal to zero also for 
H ^ 0. In the magnetic field directed along the z-axis, 
the resonance frequencies are determined by the formula 
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FIG. 10: (color online) The calculated ESR spectra of the 
planar spiral antiferromagnet with uniaxial anisotropy and 
parameters of LiCu2 02 (cjio = 320 GHz, UJ20 = GHz and 
^ — = 1.4). Solid lines show the spectra. Top panel: 



H||a. Bottom panel: H ± a. Dotted line indicate the para- 
magnetic resonance spectra. Insert to the bottom panel: the 
scheme of oscillations of the spin plane at H||c. 



The reorientation transitions occurs in the following 
fields: 

If X-L > X|| and H || z 



Hrz — 



mm(wio,W2o) I xT 



7 V X± - X 

And if x± < X|| and H directed along x or y 
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respectively. In the first case, the oscillation frequencies 
in fields above the field of the spin-fiop arc given by 
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dcpending on whether the constant B is positive (n |j x at 
H > Hcz) or negative (n || y in iJ > Hcz)- In the second 
case, the frequencies are determined by expression (1), 
into which is necessary to substitute 
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osciUations arc found by numerical calculations of the 
equilibrium orientation of the order parameter, relative 
to the crystal axes and the frequencies of small oscilla- 
tions. 

Fig. [TU] shows the calculated spectra for LiCu202 in 
the model of planar spiral antiferromagnet with uniaxial 
anisotropy and parameters: ojiq — 320 GHz, a;2o = 
GHz and , /-^i- = 1.4. 



for H II x, H > Hrr and 
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A-B 

X± 



cy 



for H II y, H > H, 

In a more complex case, where not all the spirals 
twisted to one side, in the expansion of Ua can also add 
member Cuy. In this case, as well as the fields are not 
directed along the axes of symmetry of the crystal, ex- 
pressions for the frequencies are too bulky and the orien- 
tation of the order parameter and the frequency of small 
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